Structure Formation in Non-singular Higher Curvature Cosmology 
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We propose novel structure formation scenarios based on a non-singular higher curvature cosmo- 
logical model. The model is motivated by the R^ coupling of a scalar field appearing in the string 
theory, and in our scenarios the universe has no beginning and no end. We give two examples 
with explicit parameter values which are consistent with present observations of the cosmological 
structure. In the first example, the origin of structures are generated as the adiabatic perturbation 
during the chaotic inflation, while in the second the isocurvature non-Gaussian perturbation from 
the superinflating era is responsible for the structure. In the second case it is possible to generate 
primordial supermassive blackholes whose scale is comparable to what is expected in the galactic 
nuclei. 
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I. INTRODUCTION 

Inflation is unarguably one of the essential ingredients 
of the modern cosmology. It gives a natural explanation 
to the flatness and homogeneity of the present universe, 
and the quantum fluctuation generated during the expo- 
nential expansion era explains the origin of the large scale 
structures which we observe today. The period of infla- 
tion, however, cannot continue from indefinite past ||i||, so 
the problem of the initial singularity remains unsolved. 

One possible mechanism of resolving the initial sin- 
gularity is the quantum tunneling scenario g, which is 
well-motivated by the quantum mechanical picture of the 
early universe. There are, however, other possibilities 
of constructing non-singular cosmological solutions, ob- 
tained by introducing higher curvature terms which are 
ubiquitous in quantum gravity and string theory. Con- 
struction of non-singular cosmological models by resort- 
ing to generalized gravity theories deserves serious spec- 
ulations since quantum gravitational effects are expected 
to dominate in the very early universe. 

In this paper we shall concentrate on the higher cur- 
vature resolution of the initial singularity, and propose 
a non-singular cosmological scenario consistent with the 
current observations. If the initial singularity were to be 
removed, the curvature appearing in any effective action 
must stay finite in the course of the cosmological evolu- 
tions [pi. In our model the initial singularity is avoided 
by the dominance of the Gauss-Bonnet term coupled to 
a scalar field. A merit of this model is the absence of 
redundant degrees of freedom which may appear by in- 
troducing higher curvature terms. The exact form of the 
B? coupling has been calculated in the context of string 
theory with certain compactification to four-dimensions 
, and non-singular cosmological solutions based on this 
gravity theory are known. The basic picture we propose 
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in this paper is that the avoidance of the initial singu- 
larity takes place before the ordinary scenario of chaotic 
inflation. We estimate the expected perturbation spectra 
in this scenario and discuss its implication to the cur- 
rent observations. We show that there are three different 
quantum fluctuation sources in our scenario: the isocur- 
vature fluctuations with parameter-dependent spectral 
index, the adiabatic fluctuation of blue spectrum with 
index n^ = 10/3, and the adiabatic perturbation of flat 
spectrum from the chaotic inflation. We shall see that the 
singularity avoidance in this model may leave its trace as 
a perturbation of order unity in the observationally in- 
teresting scale, which may end up now as the black holes 
at the center of the galaxies. 

The rest of this paper is organized as follows: We 
briefly review the non-singular cosmological model pro- 
posed by Antoniadis, et.al. g] in Sec. II. In Sec. Ill 
we present the outline of our non-singular cosmological 
scenario. The perturbation spectra in this model is dis- 
cussed in Sec. IV, and the structure formation scenarios 
with explicit parameter values are discussed in Sec. V. 
In Sec. VI we conclude with some comments. 



II. RESOLUTION OF INITIAL SINGULARITY 
BY HIGHER CURVATURE 

We start with an action of higher curvature gravity 
motivated by a string loop correction, 




S = 



where R%b 



i(v^)' 



-^S,{'^)RqB + ^matter 
(1) 



p JDflb'KX 



ARa^R^"' + i?2. We assume 



the function ^((/j) to be the one appearing in the string 
loop corrections H], 



^(ifi) = -ln[2eV(Je'^)], 



(2) 



which is an even function of ip, and t] is the Dedekind 
r] function. The field (p is the modulus appearing along 
with the orbifold compactification. We neglect the contri- 
bution Cmatter froni Other matter contents for a moment. 
The cosmological solutions from this action is studied in 
Variation of the action (m using the flat Friedmann- 



Robertson- Walker metric ds^ = —dt^ 



Sijdx'^dx^ , 



= 0, 
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0, 


(5) 



where a is the scale factor, yields the equations of motion 

{2H + 5H^){1 - Hi) + H^{1 
(p + iHif + 'i{H + H'^)H^tv 

where H denotes the Hubble parameter. It can be shown 
that the system is invariant under ip <-> —if, and that the 
signs of H and Lp are conserved. In the following we 
assume for definiteness the initial ip to be negative and 
'P>Q. 
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FIG. 1. The cosmological solution obtained from the action 
(qI). The universe makes a smooth transition from accelerat- 
ing phase to Friedmann-like phase. After the transition the 
scalar field loses its kinetic energy and slows down. 

Cosmologically interesting solution which leads to a 
Friedmann-like universe in late time is such that H > 0. 
This solution describes the universe which starts from a 
Minkowski space with large enough \ip\ at i = — oo, grad- 
ually accelerates its expansion and increases the moduli 
kinetic energy as 






H ^ T-TTT, if ^ ipp — 5hl\t 



^^vr' 



(6) 



where ui, ap, and ipp are constants. The constraint equa- 
tion (P) gives a relation be'^p = 2a;f tt between the initial 



values, so the initial degree of the freedom is only one, 
apart from the overall scaling of Up. Note that, reflect- 
ing the special choice of the higher derivative terms (the 
Gauss-Bonnet combination), there are no redundant de- 
grees of freedom which may make the phase space un- 
physically large. Subsequently, the Hubble parameter 
reaches its maximal value H = H-a^i^y^ at (^ ~ 0, where 
moduli loses its kinetic energy and slows down. Af- 
ter the Hubble peak, the universe naturally goes into a 
Friedmann-like phase, where 
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(7) 



Here, a/ is a constant. An example of the numerical solu- 
tion is shown in Fig. 1. In this model, there is no horizon 
problem since there is no particle horizon. It is worth 
mentioning that this model is not plagued by a run- 
ning scalar field, since the scalar field slows down in the 
Friedmann-like phase, which is similar to the Damour- 
Polyakov mechanism [|6| . 



III. COSMOLOGICAL EVOLUTION 

Realistic cosmological models need to explain the ther- 
mal equilibrium in the early universe and the generation 
of present cosmological structures. Following the usual 
inflationary scenario, we assume that the reheating is 
caused by the coherent oscillation of the inflaton at the 
end of the chaotic inflation which started shortly after 
the Hubble peak in the model described in the preced- 
ing section. This chaotic inflation is caused by the large 
fluctuation around the Hubble peak which kicks up the 
inflaton to sufficiently large expectation value. 

As a component of the matter Lagrangian, we assume 
the infiaton Lagrangian of the form, 
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inflaton 
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where for simplicity we only consider the potential 



(8) 
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After the Hubble parameter reaches its peak, the inflaton 
expectation value becomes large due to the large fluctu- 
ation. The universe region with sufficiently large x satis- 
fies the slow-roll condition and inflates. A natural initial 
condition for the chaotic inflation is estimated as 
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The subscript / denotes the value at the beginning of the 
chaotic inflation, which is supposed to be shortly after the 
Hubble peak. Then the e-folding number of the chaotic 



inflation is determined by the initial value of the inflaton 
as 
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My 
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In our scenario there are four distinct eras. The first is 
the superinflationary era, which we call SI for short, with 
increasing Hubble parameter and accelerating modulus 
field. The second is the chaotic inflationary era, denoted 
CI, where the Hubble parameter stays constant and the 
scale factor increases exponentially. At the end of this 
chaotic inflation the energy of the oscillating inflaton is 
converted into radiation, and the universe becomes ra- 
diation dominated (RD). As the radiation is redshifted 
away more quickly than the matter, in the course of time 
the universe will be matter dominated (MD). 

IV. PERTURBATION SPECTRA 

In our scenario, the inhomogeneity of the universe, 
which is expressed as the perturbation on the homoge- 
neous universe, is generated as the quantum fluctuation 
inside the horizon. The quantum fluctuation is once red- 
shifted out of the horizon and classicalized in SI and CI 
eras, and reenters inside the horizon in RD and MD eras. 
Since this inhomogeneity is supposed to be the seeds of 
the large scale structure of the present universe, the spec- 
trum of the fluctuation is crucial in predicting the large 
scale structure today and thus in discussing the plausi- 
bility of this scenario. In this section, we speculate on 
the quantum fluctuation spectra generated during the SI 
and CI eras. 



A. Adiabatic perturbation from superinflation 

In order to describe the fluctuations we expand the 
metric and the modulus up to first order of perturbations: 

ds^ = -(1 + 2(j))dt^ + 2aV\idtdx^ 

W {(1 - 2^)(5„ -I- 2E\,j) dx'dx', (12) 

ip = ip^"^ + Sip. (13) 

For the study of the fluctuations in SI era, it is convenient 
to use the uniform field gauge {5(f = 0, E = 0). In this 
gauge the perturbation of the Einstein equation reduces 
to one wave equation 0| 
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where 



A = 



B 



a'-^a{ba'^ - 2a -|- 1) 
(3a -1)2 ' 

(a-l)2(3r-2) 
5a2 - 2a + 1 



(14) 

(15) 
(16) 



Here, a and T are defined as a = 1 — H^ and F — 
—2H/3H^, respectively. As can be seen from (||), a is 
(proportional to) the fraction of the modulus kinetic en- 
ergy to the geometrical kinetic energy. F is the effective 
adiabatic index F — 1 = —G^i/SG'^o, where G^i, is the 
background Einstein tensor. They behave as a ~ t^^ 
F ~ i in the past asymptotic region and a ^ 1, F ~ 2 in 
the future asymptotic region (in the sense of (0)). Since 
the scale factor is almost constant in the past region, 
a^a rapidly decreases as the universe evolves towards the 
Hubble peak. This leads to a decrease of the effective vol- 
ume A and the perturbation behaves as if the universe 
is collapsing. Note that the other variables 4> and V are 
written in term of V', as 
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' aa [{'ST ~ 5)^ - 4)] dt. 



(17) 
(18) 



In the ordinary inflationary universe scenario, the fluc- 
tuation is quantum mechanically generated inside the 
horizon and is classicalized as it is stretched to super- 
horizon. We assume the similar picture in SI era as well 
as in CI era. According to the conventional prescription, 
we regard tp as an Schrodinger operator and decompose 
it into mode functions 

1/;= /"d3fc(27r)-3/2/a^uj^(t)e'*=^ + atu*(i)e-^'=4, (19) 

where a^ and aj. are annihilation and creation operators 
satisfying the commutation relation [dk, dU — S'^{k — q). 

The mode functions u^ and u^ are determined so that ^ 
satisfies the equation (Ol) . Using the past asymptotic so- 
lutions (0) in the superinflationary regime, the equation 
for Uk becomes 



Uk 
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4tk^ 
5ujia] 



:Uk = 0. 



(20) 



The canonically normalized solution of this equation 
which reduces to a positive frequency mode in the short 
wavelength limit (fc ^ aH) will be 

1/3 \ 

t] +iA[ 
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^5uJia^ 




where Ai and Bi are the Airy functions and Nk 
normalization factor 



Nk 



^1/2 



-1/3 
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(22) 



The long wavelength limit (fc ^ aH) of this positive 
mode function represents the behavior of the perturba- 
tions outside the horizon. The power spectrum estimated 
shortly after the Hubble peak has the form 



P..|iter~O(10-.)x4(i)'(^) , (23) 

where the asterisk denotes the values estimated at (or 
more precisely, just after) the Hubble pealqj. This 
is a steep blue spectrum with a spectral index n^ — 
dlnP^/dlnk + 1 = 10/3. It should be noted that this 
expression includes no independent variables, and is de- 
termined only by one background initial condition lui. 
Thus the amplitude of the fluctuation is determined un- 
ambiguously once the background is fixed. 



B. Isocurvature perturbation from superinflation 

Apart from the adiabatic fluctuation described above, 
there is another possible source of the fluctuation with 
significantly different spectrum in SI era. We consider a 
matter Lagrangian of an axionic field |9| , 



f^raatter — ^^ 



2tiip 



[Vaf 



(24) 



We decompose the axion field into its background part 
and a small perturbation, a — ct*-°^ -(- 5a, and suppose 
CT*^°^ = so that the background cosmological solution 
and the perturbation ip are unaffected by the introduc- 
tion of this axionic field. The equation for 5a is 



C. Perturbation from chaotic inflation 

Finally, there is a perturbation of flat n — 1 spectrum 
arising from the chaotic inflation. The amplitude of the 
density fluctuation at the horizon is estimated as 



chS^ 



2j|gl/2^3/2^ 



(27) 



Here, Nx is the e-folding number of the chaotic inflation 
since the scale A under consideration exits the horizon, 
until the time of reheating. Similarly, C,x is the value for 
the corresponding comoving fluctuation scale A. 



D. Fluctuation spectrum in our scenario 

To sum up what have been calculated, we have three 
sources of perturbations with different spectra: 

1) Adiabatic perturbation of n^ — 10/3 blue spectrum 
from superinflationary era, 

2) Isocurvature perturbation of nsa = 3 — 10/^ spectrum 
from superinflationary era, 

3) Adiabatic perturbation of n^^ = 1 flat spectrum from 
chaotic inflationary era. 

InP 
k 



5a + {3H + 2fi(p)5a 



-5a ^ 0, 



(25) 



and the canonical quantization of 5a using the past 
asymptotic background solution (o) yields the spectrurm 
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2i°^+ir(5M + ^)'. 



(26) 



The spectral index, ns^ = din Pg^j/ dink +1 = 3 — 10/i, 
depends on the coupling /i. For fi — 0.2, this gives a scale 
invariant flat spectrum. It is noted that the statistics of 
this fluctuation is non-Gaussian. 



* The adiabatic perturbation ip in the large scale limit does 
not stay constant in this model, but </!> remains constant ex- 
cept in the vicinity of the Hubble peak. We estimated the 
power spectrum (£31) by using the constancy of on the back- 
ground (0) and (IWj and we discarded the constant mode of 
the long wavelength limit solution. The classical evolution of 
perturbations is discussed in S. 

t Unlike ^, Sa stays constant outside the Hubble horizon as 
long as /i > 0. 
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FIG. 2. Perturbation spectra in our scenario. Sa is the 
isocurvature perturbation and tp is the adiabatic perturba- 
tion generated in the superinflationary era. In this figure the 
perturbation Sa is drawn as flat, but this can be tilted to 
red or blue, depending on the axion type coupling constant 
fi which is chosen accordingly. The wavenumber correspond- 
ing to the horizon scale at the Hubble peak is denoted by fc* . 
The other indicated wavenumber ko corresponds to the scale 
where the amplitude of ip becomes comparable to that of Sa. 

There might exist other sources of perturbations, such 
as the isocurvature mode perturbation during CI era. 
However, we neglect them and assume these three pertur- 
bations dominate for different scales and determine the 
whole primordial density perturbations. Since the uni- 
verse is always expanding, the perturbation generated in 
earlier time is expanded to larger scale at present. Fig. 
shows the expected perturbation spectra in our model. 



where Scr and ij) is the isocurvature and adiabatic per- 
turbations, respectively, in SI era and 5x is the adiabatic 
perturbation in CI era. The inchnation of 5a can be tilted 
in either blue or red direction depending on /i. The figure 
shows the flat case where fi = 0.2. 



V. STRUCTURE FORMATION HISTORY 

Let us discuss the observational viability of our sce- 
nario. In this model there are two distinct scales in fluc- 
tuations. One is the scale corresponding to the hori- 
zon scale of the Hubble peak, which we denote by A*. 
The comoving wave number corresponding to this scale, 
fc» ~ 27ra/A,, is shown in Fig. || The other scale, Ao, is 
where the amplitude of the blue adiabatic perturbation 
(-0) is comparable to that of the isocurvature perturba- 
tion {5a). The corresponding comoving wave number 
ko — 27ra/Ao is also shown in Fig. g. These scales are de- 
termined by the reheating temperature Trh and the initial 
inflaton value x/: through the e- folding number during 
and after the chaotic inflation up to present. Here we 
discuss two examples with explicit parameter values. 

Scale 
i 




Time 



FIG. 3. A schematic picture of the Hubble horizon scale 
Ih and the perturbation scales in our scenario. Shortly af- 
ter the Hubble peak the inflaton acquires a large expectation 
value owing to the large fluctuations, and the chaotic infla- 
tion starts. At the end of the chaotic inflation the universe 
is heated by the coherent oscillation of the inflaton. See Sec. 
VB for the explanations of the scales Igai and Ipbh- 



A. CI structure formation scenario 

We first consider the case where the scale A* is large 
enough compared to the horizon scale today, so that all 
of the observable structures are generated during the 
CI era. In this case our scenario is nothing but that 
of the ordinary chaotic inflation, though the history be- 
fore the onset of the chaotic inflation is included in our 
picture. Since the evidence of the singularity avoidance 
(the physics around the Hubble peak) has been inflated 



away beyond the present horizon, this model is consis- 
tent with our current observations only if the chaotic in- 
flation successfully takes place and generates appropriate 
amount of fluctuations. In order to predict the fluctua- 
tion amphtude of 0(10"^), the couphng constant has to 
be g ~ 0(10^^^). If we assume the constant Hubble pa- 
rameter during the CI to be 10"'' Mp, then the condition 
@ imphes xi - 4.9Mp. From (|n|), the e-folding num- 
ber of the CI era will be iV ~ 74 and we can see that 
the present horizon scale ('^ Gpc) is generated within 
the CI era. Thus, with appropriate choice of parameters, 
the chaotic inflation naturally takes place and generates 
fluctuations with observationally acceptable amplitude. 



B. SI structure formation scenario 

As another possible scenario of structure formation, 
we can consider the case where both Ao and A* are inside 
the horizon today. In this case, what generates the obser- 
vationally supported nearly flat spectrum which evolves 
to the present large scale structure, is the isocurvature 
perturbation pQ] of SI era with appropriate choice of 
/i . Then the outstanding pointed spectrum of the ip 
perturbation and the flat spectrum of 5x perturbation, 
as is shown in Fig. 0, come inside the present horizon. 
This may fall into the scale of observational interest. We 
choose, for example, the parameters as 



Trh 



O.lAfp, Ai = 0.12, g^ 
= 10^2 Qgv, and Hci 



10 



-12 



XI = 4.6Mp, 



Hi = IQ-^Mp. 



(28) 



Here, Hci is the energy scale of the chaotic inflation, 
which follows from the condition (EG). Using (O), the 
e-folding number of the chaotic inflation is estimated as 
N ~ 67. Taking account of the reheating temperature of 
10^^ GeV and rather high energy scale 10~^Mp during CI 
era, the large scale fluctuation of the universe which we 
observe today as CMB corresponds to the scale which ex- 
its the horizon during SI era. From fi — 0.12, the isocur- 
vature perturbation during the superinflationary era has 
a blue tilted spectrum with ng^ = 1.8. The amplitude of 
the fluctuation for the present horizon scale is, from (Eq) 



1/2 



10- 



with LUi ~ 0.1 and a^/ap ~ 5, estimated as Pg^ 

The fluctuation generated during the chaotic inflation 
comes to scales smaller than kpc today. If the amplitude 
of this small scale perturbation is too large, it will pro- 
duce primordial blackholes. The primordial blackholes 
of masses smaller than 10^^ g will Hawking-evaporate 
within the age of the universe, and if the radiation emit- 
ted through the evaporation is too strong, it will spoil 
the observationally supported nucleosynthesis scenario. 
Thus the amplitude of this small scale perturbation is 
constrained by the amount of the primordial blackholes. 
Assuming the spherical collapse model with Gaussian 
perturbation, the primordial blackhole mass fraction at 



the time of their formation is estimated in terms of the 
perturbation ampUtude d{M) for the mass scale M as 



f3{M) ~ S{M) exp 



1 



18(52 (Af) 



(29) 



The strongest constraint for the amount of the smaU 
blackholes is P/{1 - /3)<10-26 ||i|. Using our choice of 
the parameters (p8|), the amphtude of the perturbation 
generated during CI era, from (|27|), is ^ ~ 10"'*. Substi- 
tuting this into ( p9| ) , we can see that the smaU primordial 
blackholes are not overproduced. 

The adiabatic perturbation of nw, = 10/3, whose maxi- 
mal amplitude is determined by (^), will produce certain 
amount of primordial blackholes on its horizon reenter. 
The mass of these blackholes is determined, according 
to the standard scenario of primordial blackhole forma- 
tion, by the Hubble horizon scale at the time of reenter. 
In our scenario the sharp peak of the fluctuation spec- 
trum exits the horizon at i* , inflated by the chaotic infla- 
tion and reenters the horizon during the radiation dom- 
inated era of the universe. For our parameters (pq), the 
scale of the perturbation peak enters the horizon where 
H ~ 10~2sec~^. Then the mass of the typical primor- 
dial blackholes generated by the n^ = 10/3 spectrum 
is M ~ IO^Mq. The amount of these blackholes can 
be estimated by using (p9|). We have i5 ~ PJ ~ 0.05 
for our parameters, so /? ~ 10"^*^ at the time of horizon 
reenter. Translating this into the present fraction of the 
blackholes with respect to the whole energy density of 
the universe, we have flsH ~ 10"^. Since the spectrum 
of the adiabatic perturbation in SI era is steeply blue, 
its contribution to the perturbation amplitude in larger 
scales is negligible. 

Recent observations suggest the existence of supermas- 
sive blackholes of masses 10^ ^ IO^Mq at the galac- 
tic center. The amount of these blackhole is estimated 
using the radiation from the quasar, which is emitted 
along with the accretion of matter onto the blackholes 
p3| . The mass density of dead quasars is estimated 
as p, = u/ec^ ~ 2.2 x 105(0.1/e)Afo/Mpc^, which fol- 
lows from the observed radiation density of quasars, 
u ~ 1.3 X 10~*^erg/cm^. In the expression of p,, e is the 
efficiency of the radiation process. Assuming that the 
mass density of dead quasars are comparable to that of 
blackholes, this indicates i^BH ~ 10~^. Thus, our model 
with parameter values (Eq) predicts both mass and den- 
sity of the blackholes which are typically found in galactic 
centers. 



model, and discussed two examples with explicit param- 
eter values. Our rough estimates show that, with ap- 
propriate choice of parameters, this model can produce 
observationally acceptable amount of perturbation am- 
plitude in large scales, and even primordial blackholes in 
interesting mass scales. 

We have just assumed that the chaotic inflation, due to 
the large fluctuation, takes place after the Hubble peak, 
and we have not discussed the physics of the onset of the 
chaotic inflation in detail. Of course there is a room 
to discuss other possibilities by including other fields, 
etc., but this is beyond the scope of this paper. Here, 
we comment on one possible alternative scenario based 
on the study of the perturbation behavior. Our back- 
ground cosmological solution depends on the initial con- 
dition with only one degree of freedom. In stead of wi 
in the discussion above, here we specify the background 
by H^, the maximum Hubble parameter. For universes 
with small H^,, B of ( [l6| ) remains negative throughout 
the evolution of the universe. For iJ, larger than unity 
in our unit, there appears shortly after the Hubble peak 
a period where B becomes positive so the perturbation 
becomes unstable in short length P]. Since our unit is 
roughly Planckian, the Hubble parameter of the thresh- 
old of the instability is Planck scale. If this instability is 
physical, the perturbation inside the horizon will grow ex- 
ponentially in time and the region inside the horizon will 
collapse to form Planckian scale primordial black holes. 
Since small scale blackholes evaporate almost instantly, 
the universe will be heated by the Hawking evaporation 
of these Planckian scale blackholes. Assuming that most 
of the energy density is converted to radiation during 
the unstable period, the reheating temperature of this 
scenario would be estimated as H^ ^ ^LfP ^ N^ffT^i^, 
indicating the reheating temperature around the Planck- 
ian temperature. Here, iVg// is the number of particle 
species at the reheating epoch. If this reheating takes 
place, we have to consider a thermal inflation rather than 
the chaotic inflation in the following scenario. This infla- 
tion is necessary in order to dilute unwanted relics coming 
out of the evaporating blackholes. This scenario, how- 
ever, might be somewhat speculative, and considering 
that the action of our starting point p) is obtained by 
the loop correction in the string theory, a sensible inter- 
pretation would be that the effective theory breaks down 
for H^O{l) and that this instability is unphysical. 



VI. DISCUSSIONS 



We have described a structure formation scenario 
based on a higher curvature non-singular cosmological 




FIG. 4. The time evolution of B (^ for different back- 
ground solutions which are specified by different values of H^ . 
Note that these are calculated on the numerical solutions of 
(H-o), and no further assumptions such as the continuation to 
the CI era are taken into account. As can be seen from wM, 
positive B means the existence of the instability for the scales 
smaller than the Hubble horizon. 
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In the uniform field gauge, the equations of motion for 
the perturbations (neglecting Tj^l^^^ter) become: 



3(3a - l)H'iJj - 6H^^(I) = 2a^yip - (3a - l)H^r-V, (All) 



(3a- l)i/0 + 2ai/; = O, 

acj) - (3r - 5)ailj + -(a^aV) = 0, 
a 

^ + 3HiP + H{1 - ^H^)^ 



(A12) 
(A13) 



= {Hi^j)- + {Hi + 2H^i + 2Hi)H<j). (A14) 
These four equations reduce to the wave equation (H). 



APPENDIX A: FORMULAE OF 
COSMOLOGICAL PERTURBATIONS 



The equations of motion derived from the effective ac- 
tion (Q) can be written as the following Einstein equation: 



G^, = r„^, + T( 



GB V \ ^matter ui 



where 



G^ = i?^ - -i?(5^, 



J- (/> 



1 



^^•^ = Lp^'-if,^ 



-ip-''(p.aS^^, 






■ aP 



rjijnatter 



1 Oi--mat 



The non-vanishing components oiTcB^v are 
Tgb = — 3iJ ^, 



T, 



GB J 



H^S, + 2{H + H^)HS, 



(Al) 



Expanding Tqb^v up to first order in perturbation gives 
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